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Abstract 

Recently a novel phase transition has been observed []] in a rotating Bose-Einstein condensate 
when the rotating frequency reaches the transverse trap frequency uj± and eventually crosses it. 
We study certain aspects of this experiment in terms of the condensation of composite bosons and 
the corresponding vortex formation using a Chern-Simon Gross-Pitaevskii theory. 
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I. INTRODUCTION 



Following the observation of a regular array of vortices in a rotating Bose-Einstein con 



densates _ (BE C) [2 1 there have been a great amount of theoretical work 



y, 10, 13,ll4,|l3 



17[ to un- 



and also some further experimental development jl6, 
derstand the behaviour of the condensate when it is rotated even faster. Apart from the 
fate of the vortex lattice submitted to a faster rotation, the interest stems also from the 
possibility of observing Laughlin-liquid like states in the condensate when the rotational 
frequency becomes almost equal to the transverse trap frequency (f2 ^ uj±)- Since the num- 
ber of atoms in the rotating condensates is very large, it has been pointed out 3, IS, I3| 
that one can still use the Thomas-Fermi (TF) approximation where the kinetic energy is 
neglected compared to the interaction energy, even though the system is very close to the 
Landau level limit. Within the TF approximation it has been shown that the vortex 

lattice may survive when Q is very close to uj±. Ho ^ reached the same conclusion using 
the Lowest-Landau level (LLL) approximation where the kinetic energy term is also frozen 
by the LLL constraint and thereby it is equivalent in this respect to the TF approximation. 
In a subsequent development Sinova et. al. analyzed the Bogoliubov spectrum in the 
LLL limit and reached the conclusion that quantum fluctuations melt the vortex lattice. 
The same conclusion was reached by Wilkin et. al. j4|l using exact particle diagonalization 
technique in a toroidal geometry. Fishcher et. al. 7| have addressed the same problem 
using an effective free energy for the interacting vortices in a rotating frame and provided a 
phase diagram consisting of regions of Laughlin liquid of bosons, Laughlin liquid of vortices 
and vortex lattice. 

In a very recent experiment at ENS, Bretin et. a/.Q reached the Landau level (LL) 
limit of a BEC. They superimpose an additional quartic confining potential above the usual 
quadratic confining potential in the transverse plane, thereby eliminating the possibility of 
centrifugal instability in the limit VL = uj^. This experiment therefore reaches a regime where 
the TF approximation does not hold anymore. The vortex lattice is stable even at a very 
high value of the rotational frequency Vt (almost 98 percent of cu^). As Vt is further increased 
and eventually taken just beyond cu^, one obtains on both sides of = uj± a blurred image 
of the vortex lattice indicating a phase transition. As VL is further increased beyond 
a condensate like region reappears in a central region of the blurred image and again new 
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vortices start nucleating. At slightly higher one recovers a collection of few such vortices 
again arranged in a regular array. 

It should be pointed out that in most cases the rotational frequency Q is equal to the 
stirring frequency QsUr- However, in the ENS experiment ^IsUr is not always equal to Q. 
Though in most cases the two are almost equal, it has been pointed out that beyond Q = uj± 
it is not correct anymore. Particularly at the highest reached value of QsUr at which a vortex 
lattice becomes clearly visible again, the frequency Q is predicted to be much lower than 
Qstir- However, it has been clearly demonstrated that when the rotational frequency of the 
condensate equals or exceeds uj±, the number of vortices in the condensate is much lower 
than expected at such a high Q. 

In the present paper we argue that those features are a manifestation of boson-vortex 



dua^ 
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ity which 
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las been earlier predicted in the framework of Chern-Simon (CS) theory 
3, Isii Is^ in order to explain certain aspects of the quantum Hall effect 
(QHE). That such a composite boson can be formed at Q = uj± has been first pointed out 
by Wilkin et.al. 0]. In the quantum Hall system, however, this boson- vortex composites 
known as composite bosons remains a useful theoretical tool to understand the existence of 
the off-diagonal long-range order (ODLRO), a property of the one body density-matrix of 
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the composite bosons rather than the one body density matrix of the electrons 
However we shall argue that for a rapidly rotating BEC, composite bosons can actually be 
formed since the vortices are the zeroes of the condensate wave function and not a ficticious 
statistical flux. We propose to interpret the ENS experiment as a signature of the formation 
of composite-boson. We note that recently a CS field theory has been proposed to describe 
rotating bosons near a Feshbach resonance where we also need to consider the formation of 
molecules asaresult of the st..o„gi„te.act,o„Q. 

This paper is organized as follows. We shall start with the energy functional of the 
vortex lattice in the rotating frame as given by Fetter Q| and show why it is energetically 
favourable to form at the frequency Q = uj± a condensate of composite bosons of almost 
uniform density in a region where the quartic confinement can be considered to be smaller 
than the other relevant energy scales. We then write the free energy in terms of a Chern- 
Simon Gross-Pitaevskii (CSGP) theory which is very similar to the Chern-Simon Landau- 



Ginzburg (CSLG) theory of the quantum Hall effect 



25[ . We shall show how this expression 



of the free energy can explain some of the new features observed in the ENS experiment 
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We shall also point out that at the point Q = u± the model is self-dual. Then we make 
some suggestions on possible experiments to probe the nature of this new superfluid phase. 
Finally we summarize our observations. 



II. GROSS-PITAEVSKII ENERGY FUNCTIONAL 

The ENS experiment is performed with = 35.10^ ^''/26 atoms in a cylindrical trap 
with uj± = ujx = ujy = 27r65.6Hz,uJz = 2'kIIHz. In addition there is a confinement potential 
of the form kr'^ with k = 6.54510~^^5'./. Q|. For convenience, we shall first omit the quartic 
term assuming that it does not play any other role than making the condensate stable. 
Nevertheless we must point out that this quartic potential on its own can add interesting 
features as pointed out by Lundh (see also 21 1) and besides it will influence the collective 



spectrum around the mean-field solution we shall provide here. 

We start with the Gross-Pitaevskii free energy of the bosons of mass M in a rotating 
frame. For the order parameter \1/ 

n2 



It has been pointed out 



+ ^Mw^z^l^P + i^/l^l^ - ^*ri.r xp^] (1) 
19 1 that the vortex lattice experiences a solid body rotation [2^ 



with a solid-body velocity Vsb given by 



Vxvsb = 2n (2) 

The existence of this solid body rotation then ensures that neither the phase hS or the 
superfluid velocity Vg = ^gVS* can be spatially periodic. The system thus behaves locally 
like a superfluid in each unit cell of the lattice, but globally as a rigid body provided the 
core size of each individual vortex is less than a lattice constant. 

Under these assumption the Gross-Pitaevskii energy functional can be written as 



F = j dV[-{vg - v^bfm' + ^(V^)' + vdri)m 



Ken*|*P + ^M..,V|vl>|2 + l<^|vl>|4] (3) 
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where Vcentir^) = -^MQ^rl. 

The Landau level limit is defined by the condition V_i{r±) = —Vcent{i^±)- 
The minimum energy configuration, as we shall show using a singular gauge transforma- 
tion leading to a composites of bosons and vortices, is achieved in the mean-field approx- 
imation by considering Vs = Vsb and |\E'| = constant, namely when the kinetic energy of 
the condensate in the presence of a vortex lattice is equal to zero. In this limit the solid 
body rotation becomes identical to the superfluid motion and we recover the superfluid order 
parameter throughout the bulk. 

This new phase of composite of bosons and vortices shows superfiuidity, namely off- 
diagonal long range order (ODLRO) [2^ for the density matrix of the composite-bosons. 
It differs from the ODLRO of the bare atomic condensate with vortices obtained at low 
rotating frequencies. We assume that the energy to bind bosons to vortices results from 
the repulsive interacting energy of the original bosons, thus leading to a renormalization 
of g. This behaviour is very naturally described through a Chern-Simon Gross-Pitaevskii 
approach which has the advantage to be independent of the explicit LLL constraint. For 
further discussion we also neglect the effect of confinement in the z direction, assuming that 
the profile of the condensate in this direction is the same as for the non-rotating case. 



III. CHERN-SIMON GROSS-PITAEVSKII THEORY 

Chern-Simon (CS) field theories have been originally proposed in order to describe the 
behaviour of charged planar matter interacting with photons whose dynamics is governed 
not only by the usual Maxwell density —^F^'^F^^ but also by the Chern Simon (CS) term 
^j^'^ F^j^a^ which gives rise to topologically massive (2 -|- l)-dimensional electrodynamics 
|22| . An important consequence of the CS term is that any charged excitation also carries a 
magnetic fiux proportional to this charge. For the low-energy, long-wavelength physics in {2+ 
l)-dimensional electrodynamics the CS Lagrangian is more important than the conventional 
Maxwell Lagrangian since it contains one less derivative. For a rapidly rotating BEC, the 
density of vortices is high and is assumed to form a vortex liquid, so that we can use the CS 
Lagrangian in order to couple it to the bosonic matter density \^\'^. 

This coupling is achieved through the Chern-Simon transformation between the boson 
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order parameter \1/ and the order parameter $c6 describing the composite bosons, 



(4) 



where 6ij is the angle describing the relative position of the i-th and the j-th bosons. Here 
m is always an even integer. In two spatial dimensions, this singular gauge transformation 
attaches an integer number m of angular momentum to each boson so that 



M 



Air, 



where 



ai 



)) 



cb 



hm 



(5) 



(6) 



and A = ri X r. Here indices refer to bosons and hence repeated indices does not mean 
summation. In the CSLG description of QHE the flux of a is known as the statistical 
flux. However here the vector potentials a are associated to the vortices already present 
in the condensate, and in the limit where the number of bosons is equal to the m times 
the number of the vortices the bosons and vortices make composites. One can physically 
view this mechanism as the re-entering of bosonic matter in the vortices which are zeroes of 
the condensate-wavefunction \1/ and thereby making the condensate irrotational again. For 
m = one recovers the original Gross-Pitaevskii description. 

The Chern Simon transformation in the second quantized language is given by 



(7) 



where the operator j{r) is given by 



J{r) = im / dr^!{r)^ Ini\og{z — z)^!{r), z = re 



ie 



(8) 



3l|. Using a modified form 



This tranformation at the mean-field level gives the phase of the Laughlin-wave function. 
Considering fluctuations around the meanfield and neglectin g hig her order terms one can 
produce only the modulus of the Laughlin wavefunction 
of the composite-boson field operator (whose auto-correlation function shows true 

long-range behaviour), Rajaraman and Sondhi j^] later generalized this transformation to 
produce the full Laughlin wavefunction at the mean-field level from a second quantized 
Chern-Simon description. The same procedure can be applied here also. 
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The Chern-Simon Gross-Pitaevskii Lagrangian of the rapidly rotating bosons is written 

as 



with 

where 6pcb = ^Ib^cb - Pa and 

Lcs = 7^ I dre^'^^a^d^a, (11) 
Zirnm J 

The term proportional to 6pcb accounts for the attractive interaction between composite 
bosons. It is mediated by the gauge-fields a and thus is long ranged unlike the interaction 
between the bosons described by the order parameter 

In the original problem of rotating bosons, the additional energy due to the applied 
rotation can be written as a sum of the self energy of each vortex and of an interaction 
term between the vortices which accounts for their short range and binary interacting po- 
tential The Chern Simon transformation maps this problem onto those of a liquid 
of interacting composite bosons where the dynamical gauge fields a play a role similar to 
the magnetic field in superconductors, by canceling the applied rotational velocity in the 
bulk of the system just like for the Meissner effect, leaving composite bosons whose in- 
teraction is described by the term proportional to 6pcb in (fTUjl where p^^ accounts for the 
fact that in the absence of applied rotation, the Lagrangian Lgp vanishes at the mean field 
level. Lcs is a Lagrange's multiplier which implements the constraint relating the rotational 
flux quanta to the original matter field density. Therefore, it does not contribute to the 
free energy which is identical to those given in when transformed back. To obtain the 
the mean-field solution we first replace the operators ^cb, ^cbi ^ Lagrangian © with 
the corresponding functions $cfe(?"), '^'*b('"), a(^) by taking its expectation value in suitable 
field theoretical state. The corresponding Euler-Lagrange equations of motion are obtained 
through the functional derivatives with respect to ^*i,{r) and the dynamical gauge fields 
a(r) (For further description we shall omit the argument r). They are 
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V X a 



h 



(12) 



2M 



1 




(13) 



(14) 



where 



n 



(15) 



3a = 



2mi 



The first of these equations expresses the boson-vortex duahty while the second is the 
Maxwell Ampere's law. The last equation describes the motion of the composite bosons 
under the combined effect of an applied rotation and of the pseudo rotation generated by 
the vortices. The classical solution of these equations which minimizes the free energy is 



This solution corresponds to a constant composite boson density and the complete can- 
cellation of the external rotation by the pseudo-rotation generated by the vortex. This is 
equivalent to the Meissner effect in type II superconductors and the field (V x a) plays 
here the role of the magnetic induction B. Such a description is possible only when the 
number of vortices is high enough so that the discrete number of singular vortices can be 
replaced by a continuous function representing vortex density. At very low temperature the 
system should be well represented by this classical solution and small fiuctuations around it. 
In a real system we expect to have density modulations and the corresponding gauge field 
fiuctuations. 

What will be the effect of an increasing rotation? When A is changed from its meanfield 
value —a, the additional rotation gives a corresponding density modulation. We shall now 
see how that will produce a vortex. The vortex solution has a density profile which is 
identical to the mean field solution at long distance. Therefore one can write 




(16) 
(17) 
(18) 



A + a 












(19) 
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In order for this solution to satisfy aymptotically ()14|). 

^ + « = (20) 

must be satisfied at r — oo. In this way the vortex solution is accompanied by a change 
in the rotational velocity due to the gauge field a. Upon integration, the right hand side of 
fl20|) gives the extra flux (-^) associated with the deviation of a from its mean-field value 
—A. This is equal to the circulation flux associated with a single vortex. Thus vortices 
in composite-bosons are like quasi-hole excitations over a Laughlin-liquid like ground state 



l25l |. The extra energy of a vortex can be evaluated by numerically solving ()14j) for a 
solution of the form /(r)e*^ which satisfies and ((201). In the CSLG theory of QHE such 
solutions have been obtained by Tafelmayer, Curnoe and Weiss In the bosons- vortices 
composites, the energy required to create a vortex is going to be different from that in a 
bare atomic condensate. This is because of the renormalisation in the interaction strength 
of the composite bosons. 

The above description does not include the effect of the confining quartic potential. It 
leads to a more pronounced effect at the edge of the system where the above description 
does not hold anymore. It will also affect the collective excitation spectrum around the 
mean-field solution. It is possible that the chiral Luttinger liquid description of QHE edge 
states may capture the phenomenology at the edge 3^ . 

IV. POSSIBLE PROBES FOR COMPOSITE BOSONS 

The Chern-Simon descriptioii provides a relation between the composite-boson density 



and the circulation flux quanta 
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28|. This explains qualitatively why we should observe 
similar features in the rotating composite-boson condensate and the rotating (bare) Bose- 
Einstein condensate. However, it has been already pointed out that vortices in composite 
bosons are statistically different from their counterpart in the boson condensate. In the first 
quantized language they are equivalent to quasihole excitations in the bosonic Laughlin- 
liquid. A Berry phase measurement like the one suggested by Paredes et. a/.[^ may reveal 
their statistics. 

A composite-boson superfiuid shows an ODLRO given in the mean field by 

G{z - z') =< $r i*i(^)*c.(^')i<fr > (21) 



which can be observed in the power law decay of the auto-correlation function of a single 
composite-boson 3^ . In the QHE this correlation function is difficult to measure since it is 
built out of a composite made of real electrons and ficticious statistical fluxes. Here however, 
it may be even possible to measure it directly as the composites are made of ordinary bosons 
and real vortices. 

There is another important difference between the ordinary BEC and the composite 
BEC. We speculate that the interaction between the composite-bosons is much weaker than 
between the original bosons. This conjecture may be verified by measuring the s-wave 
scattering length when these composite bosons are formed. This may serve as an important 
step to probe the composite boson superfiuid. Alternatively, another possible way to probe 



36| and 



the composite boson superfiuid is to look at the collective excitation spectrum 
compare it with the Bogoliubov spectrum of the bare condensate. The measurement of 
collective excitations may also reveal the renormalisation of the interaction strength between 
the composite bosons. 

For Q > uj± one can again use the Thomas- Fermi approximation for the composite boson 
condensate, but with the renormalized interaction strength. By increasing Q further than 
what has been achieved in the ENS experiment, one may be able to obtain some information 
about this renormalized interaction strength from the evolution of the optical thickness of 
the atomic cloud after a time of fiight measurement and taking the s-wave scattering length 
as a variational parameter in the resulting data In this limit, an increasing Q gives rise 
to an effective confinement potential of the form {—ar"^ + kr'^) with a > (Mexican hat) 
and as a result the size of the condensate region must grow. The ENS experiment supports 
this observation. In this region the composite-bosons are again in the lowest Landau level 
and their kinetic energy can be neglected compared to the other terms in the expression of 
the free energy. It will be interesting to study the collective excitation spectrum using the 
Thomas-Fermi approximation in this regime and to compare it to the one obtained within 
the CSGP framework 13 fil l. 

n 

At a frequency Q = uj±, the ENS experiment [1| showed that the radial distribution of the 
cloud developed a shallow local minima around the point r = where the mean-field CSGP 
approximation is better justified. By reducing the strength k of the quartic confinement, this 
region will be enlarged. In that case one shall recover superfluidity and vortex generation 
for Q > uj± over a larger area compared to the present experiment. 
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V. SUMMARY AND OUTLOOK 



We have proposed a Chern-Simon Gross-Pitaevskii theory to describe the behaviour of a 
rapidly rotating condensate at f2 > u!±. The central idea in this formulation is the duality 
between bosons and vortices. We have studied certain features of a recent experiment ^| in 
terms of this theory. The description in terms of composite bosons is valid in the limit where 
the number of vortices in the condensate is of the same order as the number of atoms. There 
is no confirmation of this aspect at least at the present stage of the experiment. However 
this requires to probe the condensate in more details around Q = uj± through an adiabatic 
changing of the rotational frequency. This should lead to a better comparison between 
the theory and the experiment. We also predict a renormalisation of the s-wave scattering 
length due to the formation of composite bosons. The statistical phase-measurement type 
experiments already suggested and the measurement of collective excitations around Q = uj± 
may give more informations about the condensate of composite-bosons, namely the ODLRO 
and the nucleation of vortices. 
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